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£N| . Abstract 

i-O . Newtonian gravity and special relativity combine to produce a gravitomag- 

rS*^ \ netic precession of an orbiting gyroscope that is one fourth as large as predicted 

by General Relativity. The geodetic effect is the same in both cases. 

CM . 

1 Introduction 

> 

The similarities between Coulomb's law and Newton's law of gravitation suggest 
that the underlying theories of electromagnetism and gravitation might be similar 
, in structure. With the completion of electromagnetic theory by Maxwell, it was 

apparent that modifications of Newtonian gravity were needed to incorporate a finite 
speed of propagation of gravitational fields. This became rather more urgent after 
the introduction of special relativity with its limitations on Newtonian mechanics. 
Although it might have been natural to seek minimalist relativistic extensions of 
Newtonian gravity, this is not the way relativistic gravity theory developed. In a 
bi), bold and ingenious stroke, Einstein placed General Relativity (GR) front and center 

as the relativistic gravity theory. 

Although it is not widely appreciated, most of the weak field tests of GR can 
be successfully passed by imposing minimalist extensions on Newtonian gravity. 
However, the Gravity Probe B experiments presently underway constitute new tests, 
(see http: / /einstein. stanford.edu for weekly progress reports.) The celebrated Lense- 
Thirring "frame dragging" of GR, which will be measured by Gravity Probe B, has 
an interesting counterpart in Newtonian gravity when it is extended for the magnetic 
effects necessitated by special relativity. Rudimentary electromagnetic theory and a 
simply extended Newtonian gravity theory are compared here and the consequences 
examined for orbiting gyroscopes. 

The following transformation equations can be found in many places, e.g., p237 
of "Special Relativity" by A. P. French [1]. They describe the transformation of 
quantities between the usual inertial frames S and S'; the latter moving at speed v 
down the common +x axis direction. 
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x' = -f(x-vt) x = j(x' + vt') t' = 7(t - vx/c 2 ) t = j(t' + vx/c 2 ) (1) 

7 = l/^/l - f 2 /c 2 y = y z = z' (2) 
Transformation of velocity components: 

u x - V U x i + V . . 

= i t~2 = 7^ — TTT ( 3 ) 

1 — vu x /c L \ + vu x /c z 

u y' = i / ? u ?/ = , , / 9 ( 4 ) 

y l-vu z /c 2 y l + vu' x /c 2 

uzh uz'h /,x 

Uz ' = i Ti Uz = 77 TT~2 5 

l — vu x /c z l + vu x /c z 

Transformation of force components: 



F ,_ F x -(v/c 2 )(F.u) F x + (v/c 2 )(F>.u>) 

l — vu x /c 2 l + vu' x /c 2 



(6) 



p = p y/^ F = F y fl (7) 

v l-vu x /c 2 y l + vu' x /c 2 {) 

p/ _ Fzh p _ F'zh ,q\ 

2 l-vu x /c 2 2 l + vu' x /c 2 {) 

To describe "magnetic" effects originating from a source of a force field, consider 
two particles. Let particle 1 be affixed to the origin of the S' system and 
let particle 2 be moving at velocity u' relative to particle 1. Let the force 
F' be the force of interaction between particles 1 and 2 and consider it as arising 
from the source, 1, and acting on 2. 

We can write out the components of F applicable in system S, starting with the 
right members of Eqs (6), (7) and (8) and then substituting from Eqs. (3), (4) and 
(5) for all components of u' to obtain 

F x = F'x+^v/c^F'yUy + F'zUz) F y = 1 F' y (l - vu x / c 2 ) F' z = 1 F' z (l -vu x /c 2 ) (9) 
Lump the quantities that are independent of u into a "static" force, F st such that 

F st = F' x \ + 7^j + jF' z k (10) 
Then define a "magnetic" field quantity, Kf , as 

K f = v x F st /c 2 (11) 
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With these definitions, it can be easily verified that the remaining velocity dependent 
terms of Eq (9) can be expressed as u x Kf . Thus we can write the force F acting 
on particle 2, as observed in S as 1 

F = F st + u x K f (12) 



At the instant t=0, let F' be a coulombs law force between charges 1 and 2. 
Then in Eq. (10) F' x = kqiq2x'/r' 3 , F' y = kq\q 2 y ! /r' 3 and F' z = kq\q2z'/r' 3 where 
k = 9 x 10 9 in MKS units. But at t=0, x' = jx. With y' = y and z' = z, all 
components of Eq (10) have a multiplier of 7 and so F st = ^kq^rjr' 3 and 

K f = v x F st /c 2 = q 2 v x E/c 2 (13) 

where 

E = F st /q 2 = jkq^/r' 3 (14) 

It should be obvious here that B = Kf/q^ is the magnetic field at r from a charge 
moving with speed v in the +x axis direction. 

This little exercise has so far been a way of getting the usual magnetic effects 
if we consider forces between moving charges, but it is not restricted to electric 
forces. We are free to consider F' to be any kind of force whatever for which we 
can consider particle 1 to be a source at rest at the origin of S'. It could be a 
Newton's law gravitational force or a Hooke's law force or a Morse potential force, 
or whatever. It also is not restricted to being a central force. A dipole-dipole force 
with orientation dependence could even be considered, though it would be extremely 
messy. The point is that "magnetic" forces are straightforwardly imposed by special 
relativity alone! 

An important application of Eqs. (13) and (14), is that we get the Biot-Savart 
law from them. All we have to do is appeal to vector superposition. Consider 
a (generally curved and closed) line of moving charge. Let dqi be an increment of 
charge on the line and d\ an increment of length along the line in the S system. Since 
charge is invariant with respect to velocity, we can write dqiv as (dq\/dt)d\ = idl 
and immediately see that 

dB = (jki/c 2 )d\ x r/r' 3 (15) 

This is the contribution to B at r from a charge moving at speed v at t = 0. Now 
it is true that if the line of current is not straight, we need multiple frames S' with 
different orientations in order to calculate each contribution dB. But this is of minor 

1 We could make another choice here to define Kf as a quantity having true units of force by 
defining Kf = v x F st /c) and then having u x Kf/c as the last member of Eq. (12). This is the 
choice made in gaussian units. 
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consequence. All of the various corresponding S frames are at rest with respect to 
each other. All we have to do here is believe that we can do a vector sum of all of 
the dB's. 2 Two other things are of interest here. First, at t = in each comoving 
S' frame, x' = 73; and y' = y and z' = z. Then we have r' = (7a;) 2 + y' 2 + z' 2 . In 
the S frame chosen for the sum of all of the vector increments, the coordinates of 
the field point , 2, are (x,y,z). Second, for ordinary currents, the drift velocity is so 
small that 7 = 1 to about 11 digits. 

2 Extensions to Gravity: 

There is absolutely nothing to keep us from using the gravitational force between 
masses mi (source) and 777,2 (test particle) to calculate a g field given by 

g = F st /m 2 = - 7 Gm ir /r /3 (16) 

And again, we take the case t = 0. We note for future reference that mi is a rest 
gravitational mass in the moving S' frame and 7712 is the gravitational mass of the 
test particle as observed in this moving frame. These distinctions hardly matter in 
the circumstances considered here as we will take 7 = 1. 

We can also define a gravitomagnetic field B g = Kf / 777-2 and have a gravitomag- 
netic field vector 

B g = v x g/c 2 (17) 
And then the complete force on 777,2 is given by 

F = 777, 2 g + 777, 2 U X B g (18) 

Here 7772 is the gravitational mass of particle 2. By considering the gravitational 
forces given by Eq (7) as observed in the S and S' frames, one finds that these force 
transformation equations require that 7772 must be the rest mass of the particle even 
though it moves relative to the S frame. 3 For the Newtonian force law to be the 
weak field, low speed limit for gravitational forces and to retain complete consistency 
with special relativity, we find that only rest mass can gravitate. The inertial mass 
is, of course, 7777. 

We also obtain the analogue of the Biot-Savart law as 

dB s = {--fGi m /c 2 )d\ x r/r' 3 (19) 

2 In the case of a straight line of current, the factors of 7 disappear from the integrated E and 
B fields. 

3 Assuming mi to be so large that it is effectively affixed to the origin of S' , consider the time 
interval required for m,2 to fall some small way down the y axis toward the origin. Integrate 
a y — udu/dy to obtain u(y) and calculate the time interval as J dy/u. The time interval will not 
transform correctly unless only rest masses are used for gravitational mass, even though both mi 
and m,2 are moving relative to S. 
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Where i m is the gravitational (rest) mass current in the moving (not necessarily 
straight) line of mass. For present purposes we can take 7 = 1 and thus we could 
set r = r' = \J x 2 + y 2 + z 2 . Again, we only need to appeal to vector superposition 
in order to apply this to the calculation of gravitomagnetic fields. 



3 Earth's Gravitomagnetic Field 

Eq. (19) can be used to calculate the gravitomagnetic field of earth at the north 
pole. Here we imagine the field point for 777,2 to be at the north pole. For simplicity, 
consider the earth to be composed of spherical shells, not necessarily of the same 
density, divide the shells into mass current rings rotating rigidly and then do the 
laborious sum of the rings' contributions to obtain 

Gin , , 

B s = --vT ( 2 °) 

where / is the moment of inertia of earth about its rotation axis, r the earth radius 
and ft the earth rotational angular velocity. 

By the analogy with electromagnetism, we expect the gravitomagnetic field to 
be dipolar. Thus we write the general earth gravitomagnetic field, according to 
Newtonian gravity and special relativity as 

_ G[3(J-f)f-J] 
B s " ^3 ( 21 ) 

where J = Ift is the angular momentum vector of earth. This reduces to Eq (20) 
at the north pole. As expected, [2] this last result is four times weaker than the 
gravitomagnetic field of GR. 

One calculation of immediate interest consists of using the gravitomagnetic field 
at the north pole of earth to calculate the gravitomagnetic effect on a Foucault 
pendulum placed there. The result of this is that relative to the "fixed stars" , the 
plane of the pendulum's oscillation will precess at the rate of B g /2 = GJ/2c 2 r 3 , in 
the direction of the earths rotation, where J = 5.86 x W^kgm 2 / s is the angular 
momentum of earth relative to fixed stars. (Note B g has units of reciprocal time.) 
Numerically, this calculates out to 54 milliarcsec/yr, which is, of course, one fourth 
as large as the 220 milliarcsec/yr predicted by GR [3]. The important thing about 
this result is that it is NOT ZERO. The factor of four arises from the curvature of 
space-time of GR. 



4 Gyroscope Precession 

The results of the historic Gravity Probe B mission should be published this year 
or next, yielding a first direct measurement of the gravitomagnetic effect of the 
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rotating earth on a gyroscope in a polar earth orbit. Considering only the Lense- 
Thirring effect for now and delaying the discussion of geodetic effects, the torque on 
the gyroscope, according to GR will be 

dS 1 „ . . 

r = - = --<B g >xS (22) 

where S = Ico is the spin angular momentum vector of the gyroscope and < B g > 
is the GR gravitomagnetic field averaged over many orbits. For polar orbits, this is 

4 times the average over Eq (21), or GJ/c 2 r 3 . Thus the GR precession rate is in the 
direction of earth's rotation in the amount of 

uj lt = GJ/2c 2 r 3 (23) 

For a circular polar orbit at 650 km elevation, this yields a precession rate of 0.041 
arcsec/yr, with the gyroscope angular momentum vector precessing in the direction 
of the earth's rotation. Detailed calculations of the predictions of GR, including 
corrections for the earth J2 quadrupole moment are available [4]. The J2 correc- 
tion is negligible for the Lense-Thirring effect and only about one part in 10 3 , but 
measurable, for the geodetic effect. 

The gravitomagnetic effect of Newtonian gravity and special relativity is easily 
calculated relative to nonrotating axes with origin at earth's center. In the extended 
Newtonian theory, each increment of mass dm of the gyroscope experiences a gravit- 
omagnetic force of dmU x B g , where U = v or bit + u and u = uj x r' is the velocity 
of dm relative to the spherical gyroscope's center of mass. The torque relative to 
the gyroscope center of mass, and that causes the gyroscope to precess is 

t = J dmr' x ui x r'x < B g > (24) 

which for a spherical gyroscope, reduces to Eq. (22) again. Using Eq (21) to 
calculate < B g >, This yields a precession in the direction of the earth's rotation of 

^recess = -(1/2) < B g >= GJ/8c 2 r 3 (25) 

which is exactly one fourth of the GR result. 

5 Geodetic Precession 

Although one may think of Newtonian gravity as something existing in flat space- 
time, there is a certain sense in which space must be curved. Light slows down 
in gravitational fields and meter sticks shrink and clocks run slow in gravitational 
fields. In effect, the metric of space-time is altered. Space is inherently curved in the 
sense that the ratio of a circumference of a circle centered on a gravitating mass to 
its measured radius will not be 2n because the parts of the measured radius closest 
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to the mass will be measured with a shrunken meter stick. Interestingly, the speed 
of light, measured with shrunken meter sticks and clocks running slow will yield the 
free space speed at any location and life will seem normal anywhere. 

To account for the classic tests of GR, such as the perihelion shift of planet 
Mercury, an extended Newtonian theory must encompass decreases in the speed of 
light in gravitational fields. This is all that is needed to account for the weak field 
results that have historically been attributed to GR. But this extension is non-trivial 
when coupled with special relativity, because of the alterations of our measures 
of length and time. In many theories of gravity in otherwise flat space-time (e.g 
Krogh's deBroglie wave refraction theory [5], Puthoff's polarizable vacuum theory 
[6], Krogdahl's flat spacetime theory [7], Rastalks minimally extended Newtonian 
theory [8], the 1958 Yilmaz theory [9], etc) the changes of measures of length and 
time are equivalent to having an exponential metric of the form: 

ds 2 = c 2 dt 2 e 2<t> - e' 2<t> {dx 2 + dy 2 + dz 2 ) (26) 

Here 4>(x, y, z) is the gravitational potential (an intrinsically negative quantity) di- 
vided by c 2 ,, with c Q the free space speed of light in regions free of gravitational 
fields. It can be shown that particles (including photons) subjected to only grav- 
itational forces follow geodesies in this metric, which thereby accounts for all of 
the previous weak-field tests of General Relativity. The exponential metric depends 
only on the fact that Newtonian potentials are specified only to within an arbitrary 
additive constant, as shown by Rastall [8]. In contrast, the event horizons of GR's 
Schwarzschild coordinates black holes depend explicitly on an absolute potential of 
<j) = 1/2, which is without parallel in the rest of accepted physics. The exponential 
metric can also be derived exactly from the principle of equivalence of gravitational 
fields and accelerated reference frames in special relativity. 

In parameterized post-Newtonian (PPN) approximation, the geodetic effect de- 
pends only on the diagonal terms of the metric tensor and on the parameter 71, where 

the expansion of the time-time metric coefficient yields gu = 1 + 2</> + 271 4> 2 + 

In the isotropic coordinates form of the metric, General Relativity yields 71 = 1. 
Since the expansion of gu = e 2 ^ = 1 + 2(f) + 2(f> 2 + it is apparent that 71 = 1 
is expected in "flat" spacetime theories. Thus they will predict the same geodetic 
effect as GR. Ciufolini and Wheeler [2] give the geodetic gyroscopic spin precession 
of Gravity Probe B as 

V dV 

Wgeo = "2? x dF + (1 + 7l)VxW (27) 

where dV/dt = a + V0, with a representing any non-gravitational accelerations 
and V the orbital velocity. The first term of Eq (28) is the usual Thomas precession 
of any vector taken along an accelerated (orbital) path. The combination of first and 
second terms comprises the geodetic effect. With 71 = 1 in both GR and minimally 
extended Newtonian gravity there should be agreement on the geodetic effect. The 
geodetic effect has been measured to about 1% accuracy prior to Gravity Probe B. 
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It is difficult to say what to expect from other tensor theories, such as the 1971 
Yilmaz theory [10]. The gravitomagnetic effects of GR cannot be calculated from 
the Schwarzschild metric. A Kerr metric is required since the gravitomagnetic effects 
arise from the off-diagonal elements of the metric tensor, of which there are none 
in Schwarzschild geometry. I have no idea how to calculate the analogue of the 
Kerr metric in the Yilmaz theory, but I suspect that its off diagonal elements will 
agree with the Kerr metric to first order, which is all that will be tested in the 
earth's field. So it is possible that GP-B might confirm both the Yilmaz theory 
and GR, but on the other hand, if a correct gravity theory is a vector field theory 
in flat space-time, then GP-B should confirm the simpleminded approach given 
here; i.e., results one fourth as large as those expected from GR. Lastly, there is a 
geodetic effect from the earth orbit around the sun that would produce about 19 
milliarcsec/yr precession of the Gravity Probe B gyroscopes [3]. This is in a direction 
perpendicular to the ecliptic, which puts it very nearly in the same direction as the 
GR Lense Thirring effect. Thus, the Newtonian gravity and special relativity result 
would combine to yield a precessional angular velocity component along earth's axis 
of about .01 + .019cos(23.5°) ~ .0275 arcsec/yr parallel to the earth's rotation axis. 

In the geodetic effect the Thomas Precession term contributes one fourth of 
the effect of spacetime static curvature and with opposite sign. These effects are 
essentially the outcomes of a Fermi- Walker transport of a spin vector around an 
accelerated path and in curved space. In the Lense-Thirring effect, there should 
be no corresponding combined effects as gravitational forces are supposed to be 
entirely encompassed within GR. A "Grateful Dead exclusion rule" 4 should apply 
here. Either the GR or Newtonian effects (or neither?) should be observed, but not 
a combination of both. 
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4 e.g., wear a shirt or a tie or neither but not both. 
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